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Abstract. We discuss the scenario where the gravitino is the hghtest supersymmetric particle and 
the long-lived next-to-lightest sparticle (NSP) is the neutralino or the stau, the charged partner of 
the tan lepton. In this case staus form bound states with several nuclei, affecting the cosmological 
abundances of ''Li and ^Li by enhancing nuclear rates that would otherwise be strongly suppressed. 
We consider the effects of these enhanced rates on the final abundances produced in Big-Bang 
nucleosynthesis (BBN), including injections of both electromagnetic and hadronic energy during 
and after BBN. We show that if the stau lifetime is longer than 10^ — 10* s, the abundances 
of ''Li and '^Li are far in excess of those allowed by observations. For shorter lifetimes of order 
1000 s, it appears that stau bound state effect could reduce the ^Li abundance from standard BBN 
values while at the same time enhancing the ^Li abundance, creating a region where both lithium 
abundances match their plateau values. 

PACS. 11.30.Pb Supersymmetry - 95.35. -l-d Dark matter 



1 Introduction 

The primordial Big-Bang nucleosynthesis (BBN) pre- 
dictions for the light elements abundance provide some 
of the most stringent constraints on the decays of un- 
stable massive particles during the early Universe [1- 
8]. This is because the astrophysical determinations 
of the abundances of deuterium (D) and ^He agree 
well with those predicted by homogeneous BBN cal- 
culations, and also the baryon-to-photon ratio rj = 
Ub/rij (X f2},h? needed for the success of these calcula- 
tions [9,10] agrees very well with that inferred [11] from 
observations of the power spectrum of fluctuations in 
the cosmic microwave background (CMB). However, 
it is still difhcult to reconcile the BBN predictions for 
the lithium isotope abundances with observational in- 
dications on the primordial abundances. The discov- 
ery of the "Spite" plateau [12], which demonstrates 
a near-independence of the ^Li abundance from the 
metallicity in Population-II stars, suggests a primor- 
dial abundance in the range ^Li/H ^ (1 — 2) x lO^^*' 
[13], whereas standard BBN with the CMB value of -q 
would predict ^Li/H ~ 4 x 10"^° [9,10]. In the case of 
''Li, the data [14] lie a factor ~ 1000 above the BBN 
predictions [15]. 

The effects of hadronic injections due to late de- 
cays of the NSP during BBN have also been studied 
extensively [3,5,6,16-21]. It has recently been pointed 
out that, if it has electric charge, the NSP forms bound 
states with several nuclei [22]. Due to the large NSP 
mass (mnsp 3> Tinucicon), the Bohr radii of these bound 
states '^"^"^nucicon ~ 1 are of order the nu- 
clear size. Consequently, nuclear reactions with nuclei 
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in bound states are catalyzed, due to partial screen- 
ing of the Coulomb barrier [23,24], and due to the 
opening of virtual photon channels in radiative cap- 
ture reactions. 

Here, we present results from a new analysis [8] 
that includes the nuclear reactions induced by hadronic 
and electromagnetic showers generated by late gravi- 
tational decays of the NSP, together with the familiar 
network of nuclear reactions used to calculate the pri- 
mordial abundances of the light elements Deuterium 
(D), ^He, ''He and ^Li. In addition, we include the ef- 
fects of the bound states when the decaying particle 
is charged. We find that for lifetimes r < lO'^ — 10^ s, 
the enhanced rates of ^Li and ^Li production, exclude 
gravitino dark matter (GDM) with a stau NSP. At 
smaller lifetimes, we see that it is the ""Li destruction 
rates which are enhanced, facilitating a solution to the 
Li problems. 

2 NSP Decays and resulting showers 
during BBN 

In order to estimate the lifetime of the NSP, as well 
as the various branching ratios and the resulting EM 
and HD spectra, one must calculate the partial widths 
of the dominant relevant decay channels of the NSP. 
The decay products that yield EM energy obviously in- 
clude directly-produced photons, and also indirectly- 
produced photons, charged leptons (electrons and 
muons) which are produced via the secondary decays 
of gauge and Higgs bosons, as well as neutral pious 
(tt*^). Hadrons (nucleons and mesons such as the 

and tt^) are usually produced through the sec- 
ondary decays of gauge and Higgs bosons, as well (for 
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the mesons) as via the decays of the heavy t lepton. It 
is important to note that mesons decay before interact- 
ing with the hadronic background [3,21]. Hence they 
are irrelevant to the BBN processes and to our analy- 
sis, except via their decays into photons and charged 
leptons. Therefore, the HD injections on which we fo- 
cus our attention are those that produce nucleons, 
namely the decays via gauge and Higgs bosons and 
quark-antiquark pairs. 

For the neutralino NSP Xj include the two- 
body decay channels x ^ G Hi and x GV , where 
Hi — h,H,A and V — Z . In addition, we include 
here the dominant three-body decays x ~^ G7* — > 
Gqq, X G7* -> GW+W-, X GW+W- and 
the corresponding interference terms. For the f NSP 
case, the lighter stau is predominantly right-handed, 
its interactions with W bosons are very weak (sup- 
pressed by powers of mr) and can be ignored. The 
decay rate for the dominant two-body decay channel, 
namely f — > Gr, has been given in [25]. However, this 
decay channel does not yield any nucleons. Therefore, 
one must calculate some three-body decays of the f 
to obtain any protons or neutrons. The most relevant 
channels are f^GT*^GZT, f^Zf*^GZT, 
f^Tx*-^ GZt and f GZt [21]. 

Having calculated the partial decay widths and 
branching ratios, we employ the PYTHIA event gen- 
erator [26] to model both the EM and the HD spectra 
of the NSP decays. These spectra and the fraction of 
the energy of the decaying particle that is injected as 
EM energy are then used to calculate the light-element 
abundances, as it is described in [8]. 

3 Bound-State Effects 

One the other hand, it has recently been pointed out [22] 
that the presence of a charged particle, such as the 
stau, during BBN can alter the light-element abun- 
dances in a significant way due to the formation nega 
tively -charged staus of bound states (BS) with charged 
nuclei. The binding energies of these states are a^Zfmi/'^ 
K. SOZfAi keV, and the Bohr radii ~ (aZim^) 
1 Z^^A^^ fm. For species such as ''He, ^Li and ^Be, 
these energy and length scales are close to those of 
nuclear interactions, and it thus turns out that bound 
state formation results in catalysis of nuclear rates via 
two mechanisms. 

One immediate consequence of the bound states 
is a reduction of the Coulomb barrier for nuclear re- 
actions, due to partial screening by the stau. Since 
Coulomb repulsion dominates the charged-particle rates, 
all such rates are enhanced. An additional effect en- 
hances radiative capture channels A2{Ai,j)X by in- 
troducing photonless final states in which the stau 
carries off the reaction energy transmitted via virtual 
photon processes. In particular, the ^He(d, 7)^Li re- 
action, which is suppressed in standard BBN, is en- 
hanced by many orders of magnitude by the presence 
of the bound states. As described in [22], the virtual 
photon channel has a cross section which is enhanced 



over that of the usual radiative capture cross section 
by '^^BB'^ (aw)~", where a is the BS Bohr radius and 

uj = is the photon energy. The index n depends 
on the type of transition multipole: for (E1,E2) tran- 
sitions, n = (3,5). To account for bound state effects, 
an accurate calculation of their abundance is neces- 
sary. To do this we solve numerically the Boltzmann 
equations (13) and (14) from [23], that control these 
abundances. 



4 Results and Discussion 

Our framework is CMSSM and mSUGRA models [27], 
where the NSP could be either the lighter stau or the 
lightest neutralino. We start presenting results based 
on CMSSM models with = 0, > and tan f3 = 10, 
showing explicit element abundance contours. In Fig. 
[T^, we show the element abundances that result when 
the gravitino mass is held fixed at 7713/2 = 100 GeV in 
the presence of stau bound state effects. To the left of 
the near- vertical solid black line at mi/2 — 250 GeV, 
the gravitino is the not the LSP, and we do not con- 
sider this region here. The diagonal red dotted line 
corresponds to the boundary between a neutralino and 
stau NSP. Above the line, the neutralino is the NSP, 
and below it, the NSP is the stau. Very close to this 
boundary, there is a diagonal brown solid line. Above 
this line, the relic density of gravitinos from NSP decay 
is too high, i.e., {m3/2/'n^NSp)^NSph^ > 0.12. Thus 
we should restrict our attention to the area below this 
line. Note that we display the extensions of contours 
which originate below the line into the overdense re- 
gion, but we do not display contours that reside solely 
in the upper plane. 

The very thick green line labelled ^Li = 4.3 corre- 
sponds to the contour where ''Li/H = 4.3 x 10~^°, a 
value very close to the standard BBN result for ^Li/H. 
There are additional (unlabeled) thin green contours 
showing '^Li/H = 2 x 10^^'' (dashed). For this case 
with = 100 GeV the ^Li abundance is never suf- 
: ficiently high to match the observed ^Li plateau for 
the same parameter values where ''Li is reduced. The 
^Li/^Li ratio is shown by the solid blue contour la- 
beled ^Li/''Li = 0.15. At large TOi/2, the contour for 
^Li/'^Li = 0.01 is shown by the thin blue line. To the 
right of this contour, including the region where ^Li 
~ 2x 10~^°, the ^Li abundance is too small. Finally, we 
show the contours for D/H = 2.2 and 4.0 x 10~^ by the 
solid purple contours as labeled. The D/H = 2.2 x 10~^ 
contour is a small loop within the ^Li/^Li loop. Inside 
this loop D/H is too small. Between the two curves 
labeled 4.0, the D/H ratio is high, but not necessarily 
excessively so. As a better illustration of our results, we 
have shaded as orange (lighter) the region where the 
differences between the calculated and observed light- 
element abundances are no greater than in standard 
BBN without late particle decays, that is ^He/_D < 1, 
^Li/'^Li < 0.15, 2.2 X 10^5 < D/R < 4.0 x 10"^ and 
^Li/H < 4.3 X 10^1". Only when mi/2 > 3500 - 4000 
GeV does the D/H abundance drop back to acceptable 
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Fig. 1. The {mip,mo) planes for Aq = 0, n > and 
tan/3 = 10. In the upper (lower) panel we use m^p = 
100 GeV (jT!.3/2 = 0.2 mo). The regions to the left of the 
solid black lines are not considered, since there the gravitino 
IS not the LSP. In the orange (light) shaded regions, the dif- 
ferences between the calculated and observed light-element 
abundances are no greater than in standard BBN without 
late particle decays. In the pink (dark) shaded region in the 
lower panel the the lithium abundances match the observa- 
tional plateau values. The significance of the other contours 
is explained in the text. In these figures we have incorpo- 
rated the bound-state effects. 



levels with good abundances for ^Li, but ^Li is now too 
small to account for the plateau. Thus, for a constant 
value of TO3/2 = 100 GeV, the bound-state effects force 
one to extremely large values of mi/2 primarily due 
to the enhanced production of ^Li, as shown by the 
orange shaded region. For this value of the gravitino 
mass, there are no regions where both lithium abun- 
dances match their plateau values. 

In Fig. [TJj, we fix m3/2 = 0.2 mo and neglect the 
bound-state effects. The choices of contours are sim- 
ilar to the upper panel. The gravitino relic density 
constraint now cuts out some of the stau NSP re- 
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Fig. 2. The {mi/2,rno) plane for mSUGRA with 7713/2 = 
mo and Ao/mo = 3 — \/3 as in the simplest Polonyi su- 
perpotential. The various contours and regions are as in 
Fig.m 



gion at large mi/2 and large mo, but allows a small 
neutralino NSP region at low m-i/2- In this case the 
constraint from ■^He/D is not very strong in the stau 
NSP region and the contour is not shown. The re- 
gion where the ^Li/^Li ratio lies between 0.01 and 
0.15 now forms a band which moves from lower left 
to upper right. Thus, as one can see in the orange 
shading, there is a large region where the lithium iso- 
topic ratio can be made acceptable. However, if we 
restrict to D/H < 4.0 x 10~^, we see that this ra- 
tio is interesting only when ^Li is at or slightly below 
the standard BBN result. Once again, we see that the 
increased production of both ^Li and ^Li excludes a 
portion of the stau NSP region where mi/2 ^ 1500 
GeV for small mo. The lower bound on mx/2 increases 
with mg. In this case, not only do the bound-state ef- 
fects increase the ^Li abundance when mi/2 is small 
(i.e., at relatively long stau lifetimes), but they also de- 
crease the ^Li abundance when the lifetime of the stau 
is about 1500 s. Thus, at (mi/2,mo) ~ (3200,400), 
we find that ^Li/^Li ~ 0.04, ^Li/H ~ 1.2 x 10-^°, 
and D/H ~ 3.8 x 10^. Indeed, when mi/2 is between 
3000-4000 GeV, the bound state effects cut the ^Li 
abundance roughly in half. In the darker (pink) re- 
gion the lithium abundances match the observational 
plateau values, with the properties ®Li/^Li > 0.01 and 
0.9 X 10-1° < ^Li/H < 2.0 X 10-^°. 

Finally, we come to an example of a mSUGRA 
model in Fig. [2] Here, because of a relation between the 
bilinear and trilinear supersymmetry breaking terms: 
Bq ^ Aq — mo, tan/3 is no longer a free parame- 
ter of the theory, but instead must be calculated at 
each point of the parameter space. Here, we choose 
an example based on the Polonyi model for which 
Ao/mo = 3 — In addition, we have m3/2 — rrio- 
The upper part of the plane, we do not have GDM. 
We note that ^Li is interestingly high, between 0.01 
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and 0.15 in much of this region. Due to the bound- 
state effects both hthium isotope abundances are too 
large except in the extreme lower right corner, where 
there is a small region shaded orange. Henceforth, the 
BBN constraints and the bound state effects practi- 
cally exclude the bulk of the parameter space of this 
simple Polonyi SUGRA model. 



5 Conclusions 

We discussed the cosmological light-element abun- 
dances in the presence of the electromagnetic and 
hadronic showers due to late decays of the NSP in 
the context of the CMSSM and mSUGRA models, in- 
corporating the effects of the bound states that would 
form between a metastable stau NSP and the light nu- 
clei. Late decays of the neutralino NSP constrain sig- 
nificantly the neutralino region, since in general they 
yield large light-element abundances. The bound-state 
effects are significant in the stau NSP region, where ex- 
cessive ^Li and ^Li abundances exclude regions where 
the stau lifetime is longer than 10'^ — lO'* s. For life- 
times shorter than 1000 s, there is a possibility that 
the stau decays can reduce the ^Li abundance from the 
standard BBN value, while at the same time enhanc- 
ing the ^Li abundance, defining a region where both 
lithium abundances match their plateau values. 
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